Spurious diapycnal mixing in terrain-following coordinate models: The problem and a solution   总被引：5，自引：0，他引：5  

Patrick Marchesiello  Laurent Debreu  Xavier Couvelard 《Ocean Modelling》2009,26(3-4):156-169

In this paper, we identify a crucial numerical problem in sigma coordinate models, leading to unacceptable spurious diapycnal mixing. This error is a by-product of recent advances in numerical methods, namely the implementation of high-order diffusive advection schemes. In the case of ROMS, spurious mixing is produced by its third-order upwind advection scheme, but our analysis suggests that all diffusive advection schemes would behave similarly in all sigma models. We show that the common idea that spurious mixing decreases with resolution is generally false. In a coarse-resolution regime, spurious mixing increases as resolution is refined, and may reach its peak value when eddy-driven lateral mixing becomes explicitly resolved. At finer resolution, diffusivities are expected to decrease but with values that only become acceptable at resolutions finer than the kilometer. The solution to this problem requires a specifically designed advection scheme. We propose and validate the RSUP3 scheme, where diffusion is split from advection and is represented by a rotated biharmonic diffusion scheme with flow-dependent hyperdiffusivity satisfying the Peclet constraint. The rotated diffusion operator is designed for numerical stability, which includes improvements of linear stability limits and a clipping method adapted to the sigma-coordinate. Realistic model experiments in a southwest Pacific configuration show that RSUP3 is able to preserve low dispersion and diffusion capabilities of the original third-order upwind scheme, while preserving water mass characteristics. There are residual errors from the rotated diffusion operator, but they remain acceptable. The use of a constant diffusivity rather than the Peclet hyperdiffusivity tends to increase these residual errors which become unacceptable with Laplacian diffusion. Finally, we have left some options open concerning the use of time filters as an alternative to spatial diffusion. A temporal discretization approach to the present problem (including implicit discretization) will be reported in a following paper.  相似文献   

Impact of wildfire on source water contributions in Devil Creek,CA: evidence from end‐member mixing analysis     

Helen Y. Jung  Terri S. Hogue  Laura K. Rademacher  Tom Meixner 《水文研究》2009,23(2):183-200

A geochemical and end‐member mixing analysis (EMMA) is undertaken in Devil Canyon catchment, located in southern California, to further understanding of watershed behaviour and source water contributions after an acute and extensive wildfire. Physical and chemical transformations in post‐fire watersheds are known to increase overland flow and decrease infiltration, mainly due to formation of a hydrophobic layer at, or near, the soil surface. However, less is known about subsurface flow response in burned watersheds. The current study incorporates EMMA to evaluate and quantify source water contributions before, and after, a catchment affected by wildfires in southern California during the fall of 2003. Pre‐ and post‐fire stream water data were available at several sampling sites within the catchment, allowing the identification of contributing water sources at varying spatial scales. Proposed end‐member observations (groundwater, overland flow, shallow subsurface flow) were also collected to constrain and develop the catchment mixing model. Post‐fire source water changes are more evident in the smaller and faster responding sub‐basin (interior sampling point). Early post‐fire storm events are dominated by overland flow with no significant soil water or groundwater flow contribution. Inter‐storm streamwater in this smaller basin shows an increase in groundwater and a decrease in soil water. In the larger, baseflow‐dominated system, source water components appear less affected by fire. A slight increase in lateral flow is observed with only a slight decrease in baseflow. Changes in the post‐fire flow regimes affect nutrient loading and chemical response of the basin. Relatively rapid recovery of the chaparral ecosystem is evidenced, with active re‐growth and evapotranspiration evidenced by the fourth post‐fire rainy season. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

Determination of residence time and mixing processes of the Ubatuba,Brazil, inner shelf waters using natural Ra isotopes   总被引：1，自引：0，他引：1  

Willard S. Moore  Joselene de Oliveira 《Estuarine, Coastal and Shelf Science》2008

Coastal waters contain elevated dissolved activities of short-lived radium isotopes, ²²³Ra and ²²⁴Ra, having half-lives of 11.4 and 3.66 days, respectively. The input of these isotopes near the coast must be balanced by decay and mixing into the open ocean, where excess activities are zero. Since the decay rate is known, in the ideal case the mixing rate may be determined from the offshore distribution of these isotopes. This study found that samples collected in June 2000 followed the expected exponential decrease with distance offshore. We assign a dispersion coefficient of 28–39 m² s⁻¹ for this study. During January 2002 and November 2003, there was not a consistent decrease of activity with distance offshore. This is likely due to the ruggedness of the coastline, where many bays and small islands interrupt simple mixing patterns. To estimate exchange rates during 2002 and 2003, we used a model based on the decrease in the ²²⁴Ra/²²³Ra activity ratio (AR) with time for samples isolated from fresh inputs of Ra. This model yielded residence times of 1–2 weeks for samples collected within 20 km of the coast. We used this residence time to calculate the flux of ²²⁸Ra (half-life = 5.7 years) to the study area necessary to maintain the enrichment relative to ocean water. This enrichment is a factor of ten greater than the flux of ²²⁸Ra expected from submarine groundwater discharge (SGD) occurring within 50 m of shore.  相似文献   

Shear and Richardson number in a mode-water eddy   总被引：1，自引：0，他引：1  

Blair J.W. Greenan   《Deep Sea Research Part II: Topical Studies in Oceanography》2008,55(10-13):1161

Measurements of stratification and shear were carried out as part of the EDDIES tracer release experiment in mode-water eddy A4 during the summer of 2005. These measurements were accomplished using both shipboard instrumentation and a drifting mooring. A strong relationship between shear intensity and distance from the center of the eddy A4 was observed with the shipboard ADCP. Diapycnal diffusivity at the SF₆ tracer isopycnal prior to and during the release was estimated from the drifting mooring to be 2.9×10⁻⁶ m² s⁻¹. Diffusivity increased by an order of magnitude to 3.2×10⁻⁵ m² s⁻¹ during the period of the final tracer survey in early September, which was similar to the value estimated from the tracer analysis for the whole experiment (3.5×10⁻⁵ m² s⁻¹, [Ledwell, J.R., McGillicuddy Jr., D.J., Anderson, L.A., 2008. Nutrient flux into an intense deep chlorophyll layer in a mode-water eddy. Deep-Sea Research II, this issue [doi:10.1016/j.dsr2.2008.02.005]].  相似文献   

Experimental determination of the surface photometric contribution in the spectral reflectance deconvolution processes for a simulated martian crater-like regolithic target     

Aurélien M. Cord  Patrick C. Pinet  Yves Daydou  Serge D. Chevrel 《Icarus》2005,175(1):78-91

  相似文献   

Kinematic wave solutions for pollutant transport over an infiltrating plane with finite‐period mixing and mixing zone     

V. P. Singh 《水文研究》2002,16(12):2441-2477

Kinematic wave solutions are derived for transport of a conservative non‐point‐source pollutant during a rainfall‐runoff event over an infiltrating plane for two cases: (i) finite‐period mixing and (ii) soil‐mixing zone. Rainfall is assumed to be steady, uniform and finite in duration, and it is assumed to have zero concentration of pollutants. Infiltration is assumed constant in time and space. Prior to the start of rainfall, the pollutant is distributed uniformly over the plane. In the first case, when rainfall occurs, the mixing of pollutant in the runoff water occurs in a finite period of time. In the second case, the chemical concentration is assumed to be a linearly decreasing function of rainfall intensity and overland flow. The solute concentration and discharge are found to depend on the flow characteristics as well as the solute concentration parameters. The characteristics of solute concentration and discharge graphs seem to be similar to those reported in the literature and observed in laboratory experiments. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

江西永平铜矿区花岗岩热液蚀变与岩石成因:矿物化学、元素地球化学和Sr-Nd-Hf同位素制约   总被引：15，自引：11，他引：4  

朱碧  蒋少涌  丁昕  姜耀辉  倪培  顾连兴 《岩石学报》2008,24(8)

永平铜矿床是位于华南钱塘江-信江断裂坳陷带中的一个大型铜矿床.矿区内出露的岩浆岩主要有似斑状黑云母花岗岩体和少量石英斑岩脉,黑云母花岗岩岩体内及接触带发育斑岩型和矽卡岩型铜钼钨多金属矿化.对岩体中黑云母的电子探针成分分析表明,似斑状黑云母花岗岩中黑云母具有富镁贫铁高钛的特点,Fe/(Fe Mg)为0.21～0.33,TiO2含量为2.12%～4.28%;而石英斑岩脉中黑云母更加富镁贫铁,Fe/(Fe MS)为0.10～0.15,TiO2含量为2.43%～3.86%.黑云母的化学组成表明花岗质岩浆形成时的氧逸度很高,与长江中下游地区其他含铜花岗质岩体相吻合.永平花岗岩体遭受了较强的热液蚀变作用,部分黑云母蚀变为绿泥石.石英斑岩中绿泥石具有富铁贫镁低钛的特征(Fe/(Fe Mg)=0.72～0.76;TiO2=0.02%～0.06%),其形成温度为139～224℃.而似斑状黑云母花岗岩中绿泥石相对贫铁富镁高钛(Fe/(Fe Mg)=0.36～0.44;TiO2=0.05%～0.36%),形成温度稍高,为229～346℃,与该矿床流体包裹体研究获得的成矿温度(220～400℃)基本吻合.Sr-Nd-Hf同位素综合研究表明,永平似斑状黑云母花岗岩具有变化较大且相对较高的εHf(t)值(-0.1～-10.3)和εNd(t)值(-5.83～-7.95),反映岩石具有明显的壳幔混合成因特征.而石英斑岩的εHf(t)值(-8.4～-12.5)和εNd(t)值(-9.93～-10.2)均稍低于似斑状黑云母花岗岩,反映其形成过程中幔源物质贡献相对较小.永平岩体的地壳端元很可能就是该区中元古代基底地层重熔的产物.  相似文献   

Magma Ascent and Crustal Accretion at Ultraslow-Spreading Ridges: Constraints from Plagioclase Ultraphyric Basalts from the Arctic Mid-Ocean Ridge   总被引：2，自引：0，他引：2  

Hellevang  Bjarte; Pedersen  Rolf B. 《Journal of Petrology》2008,49(2):267-294

Plagioclase ultraphyric basalts (PUBs) with up to 54% plagioclasephenocrysts were dredged in the rift valley and adjacent flanksof the ultraslow-spreading Mohns and Knipovich ridges. The PUBsshow large variations in crystal morphologies and zoning. Thelarge variations suggest that single basalt samples containa mixture of plagioclase crystals that aggregated at differentlevels in the magma conduits. Resorbed crystals and repeatedreverse zones suggest that the magma reservoirs were replenishedand heated several times. Thin concentric zones with melt inclusions,and sharp reductions in the anorthite content of 3–7%,are common between the reverse zones. These zones, and skeletalcrystals with distinctly lower anorthite contents than massivecrystals, are interpreted to be the result of rapid crystalliztionduring strong undercooling. The changes between short periodsof cooling and longer periods with reheating are explained bymultiple advances of crystal-rich magma into cool regions followedby longer periods of gradual magma inflow and temperature increase.The porphyritic basalts are characterizd by more depleted andmore fractionated compositions than the aphyric basalts, withlower (La/Sm)_N, K₂O and Mg-numbers. This relationship, and theobservation that PUBs are sampled only close to segment centresalong these ridges, suggests that the PUBs formed by higherdegrees of melting and evolved in more long-lived magma reservoirs.We propose that the zoning patterns of plagioclase crystalsand crystal morphologies of these PUBs reflect the developmentand flow of magma through a stacked sill complex-like conduitsystem, whereas the aphyric equivalents represent later flowof magma through the conduit. The formation of voluminous higher-degreemelts may trigger the development of the magma conduits andexplain the generally depleted compositions of PUB magmas. KEY WORDS: basalt; mineral chemistry; MORB; magma mixing; magma chamber; major element  相似文献   

Different Generating Mechanisms for the Summer Surface Cold Patches in the Yellow Sea     

Minghai Huang  Hui Wu  Yihe Wang 《大气与海洋》2018,56(4):199-211

ABSTRACT

It has been frequently observed that there are surface cold patches (SCPs) in the Yellow Sea in summer. Although previous studies based on monthly mean temperature distribution found that these SCPs are a result of tidal mixing and tide-induced upwelling, tidal mixing and upwelling alone cannot explain all the occurrences. In our study we found that the three typical patches, namely, the Shandong SCP, the Subei SCP, and the Mokpo SCP, have different temporal patterns over a spring–neap tidal cycle; hence, they have different generating mechanisms. Based on a multiple-year simulation, the latter two show conspicuous spring–neap variations. The highest temperature occurs during the neap tide phase (about two days after a quarter moon). Because of weak upwelling and mixing, strong stratification is established and the SCPs are suppressed or even disappear. The opposite holds for the spring tide phase (about two days after a new or full moon). This is quite different from the Shandong SCP, which does not display a distinct difference between spring and neap tides. Buoy observations and composite analyses with data from the Moderate Resolution Imaging Spectroradiometer (MODIS) further support this conclusion.  相似文献   

基于非充分掺混模式的流域来水组成模型          下载免费PDF全文

王船海  朱琰  程文辉  向小华 《水科学进展》2008,19(1):94-98

针对如何减小数值求解对流输运方程耗散误差的问题,引进断面计算浓度的概念来计算河段平均浓度,提出了非充分掺混模式的有限控制体积法离散对流输运方程的新算法,据此构建了模拟流域内的水源组成以及不同水源在流域内的时空变化情况的流域来水组成模型。通过数值试验与具体实例验证,结果表明所提出的非充分掺混模式的新算法可有效地提高流域来水模型的计算精度。  相似文献